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ABSTRACT 

We perform two-dimensional hydrodynamic simulation on the gas stripping by radi- 
ation drag from an interstellar cloud moving in uniform radiation fields. To properly 
include relativistic radiation drag, the radiation hydrodynamic equation is solved with 
taking into account the dilution of radiation fields by optical depth of the cloud. As 
a result, it is found that the optically-thin surface layers are effectively stripped by 
radiation drag from an optically-thick gas cloud, and simultaneously stripped gas loses 
momentum. The momentum loss time-scale is found to be on the order of 10* years 
under intensive radiation fields which are expected in the early phase of galaxy evo- 
lution. The present results show that the radiation drag is an effective mechanism to 
extract angular momentum from interstellar medium and allows it to accrete onto the 
galactic center. The mass accretion driven by radiation drag may lead to the formation 
of a central supermassive black hole. 

Key words: methods: numerical - galaxies: bulges - hydrodynamics - galaxies: nuclei 
- galaxies: starburst - black hole. 



1 INTRODUCTION 

If the matter moves in radiation fields, it absorbs the radia- 
tion anisotropically in comoving frame and then reemits the 
radiation isotropically, so that the gas is subject to drag force 
in proportion to vE, where v is t he peculiar velocity and E 
is the radiation energy density jMihalas fc Mihalasl HqsI . 
This effect is known as the radiation drag. A famous exam- 
ple of radiation drag is the Foynting-Robcrtson effect on a 
dust grain in the Solar System (Povntine 1903; Robcr tsonI 
^sj), where, due to the angular momentum loss by the ra- 
diation drag, a grain spirals into the Sun. 

The importance of radiation drag has been considered 
for other various issues in astrophysics. In the epochs shortly 
after the recombination, the cosmic background radiation 
can exert strong drag force thro ugh electron scattering, 
which is called the Compton drag. iLoebl lll993l) considered 
the possibility that the Compton drag could extract angular 
momentum from an early formed object and build up a mas- 
sive b lack hole. In this context, lUmemura. Loeb fc Turiieil 
il993h performed numerical simulations, using a three- 
dimensional hydrodynamics-dark-matter code, and found 
that the angular momentum extraction actually takes place 
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efficiently at redshifts of z > 300, leading to the formation 
of massive black holes. 

The radiation drag induced mass accretion is also 
considered in a rotating gas disk immersed in the ex- 
ternal radiation fields. Cosmological accretion disks em- 
be dded in background radi at ion are extensively studied 
bvlFukue fc Umemural Jl994) . l Tsuribe. Fukue fc Umemural 
|l994l) . lUmemu ra fc F ukud (ll994^. iTsuribe fc Umemura 
Jl997li . and iMineshi ee. Tsuribe fc Umemural (1199811 . ~ 
relation to quasi-perio d ic ph e nomena in X-ray stars, 
Fortner. Lamb fc Milleij iLambl (|l98flL Il99ll ) and 



^Uer^^Lemtbri^g^ "^ examined the role of radiation 



(Compton) drag exerted by the radiation fields of a cen- 
tral neutron star on the d i sk ga s in the binary regime. 
iMever fc Mever-Hofmeisteil il994r) proposed siphon pro- 
cesses in cataclysmic variables, where the evaporating gas 
from the disk surface infall s toward a central white d warf 
because of radiation drag. iFukue fc Umemural lll995ll in- 
vestigated a steady accretion driven by the radiation of 
a central luminous object. Also, in active galactic nu- 
clei (AGNs), the radiatio n drag may play an impor- 
tant r ole in some situations. lUmemura. Fukue fc Mineshigj 
1^93) proposed the starburst-induced fueling to an AGN 
("radiative avalanche"), to relate the circumnuclear star- 
burst with AGN even t. This model has been extended to 
more realistic versions llFukue. Umemura fc Mineshig ell997l : 
lUmemura. Fukue fc Mineshigd Il998t lOhsuea et alJ Il999ll 
Commonly, it was found that the surface layer of rotating 
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gas irradiated by intensive starlight can lose the angular mo- 
mentum by radiation drag, which results in an avalanche of 
the layer as an inevitable consequence. 

The model of mass accretion induced by uniform back- 
ground radiation is revitalized in galactic bulges, because 
bulge stars produce nearly uniform radiation fields. This can 
provide a possible way to cause the formation of a supermas- 
sive black hole (SMBH) in a bulge. The recent observations 
have suggested that the mass of a SMBH tightly correlates 
with the mass of the host bulge or the velocity dispersion 
of bu lg e star s jKormen dv fc Ri chstonelll995l: iMivoshi et al " 
199,4 iLaoil Il99^ : iMagorrianI ll99Sl: iRichstone et al 
1998"; 'H d Il999l: 'Wandcl 199 3: iKormendv fc Ho 2000* 
Gebhardt et al.l .2000a b: Ferrarese fc Merritt. .200i 



Merrifield. F( 



Salucci et al. 



20011: iMerr itt & Forr 



brbes fc TerlevichI I^OOOl: iNelsonf 
."200rt 'Ferrarese et al."200l':'McL urefc Dunlor 



2001a b; 



McLure fc Dunloni l2002l : IWandell l200^ . AU of these cor- 



relations imply that the formation of a SMBH is likely to 
be physically linked to the format i on of a galactic bulge 
that harbors a SMBH. lUmemural l|200ll) considered the 
accretion of dusty gas driven via the radiation drag by 
bulge stars and const r ucted a quasar formation model. 
iKawakatu fc Umemural l^2002^ quantitatively clarified by ra- 
diative transfer calculations that the inhomogenity of dusty 
interstellar medium controls the efficiency of radiation drag, 
and found that, if surface gas layers are stripped effectively 
by the radiation drag from optically-thick dusty gas clouds 
and lose angular momentum, a sufficient amount of gas can 
accrete to form a SMBH in realistic situations. The model 
of radiation drag-induced formation of SMBH predicts the 
SMBH-to-bulge mass relation to be 0.3-0.5£, where e is 
the fraction of mc^ liberated by hydrogen burning in stars. 
However, the processes of gas stripping by the radiation 
drag from an optically-thick cloud have not been elucidated 
in detail. The study requires a sophisticated numerical 
simulation with relativistic radiation hydrodynamics. In 
this paper, we perform a radiation hydrodynamic simulation 
with generalized coordinates for a gas cloud moving in 
radiation fields. Through this simulation, we attempt to 
demonstrate the stripping efficiency and the time-scale of 
angular momentum loss. 

The paper is organized as follows. In Section |21 we 
present the model used in this work. In Section |3 we de- 
scribe the numerical treatment. In Section |31 we show the 
numerical results on the flow patterns and then estimate 
the time-scale for momentum loss of the gas. Section |S| is 
devoted to conclusions. 



2 MODEL 

We suppose a simple situation that a spherical cold gas cloud 
moves in uniform radiation fields with a constant velocity. 
The gas cloud is assumed to have a uniform density distribu- 
tion and to be in pressure balance with hot rarefied ambient 
gas. If there is relatively large velocity difference between the 
gas cloud and the ambient matter, the cloud surface layers 
are expected to be stripped by the ram pressure. The ram- 
pressure stripping itself does not extract angular momentum 
from the total system including the cloud and ambient mat- 
ter. But, if the system is embedded in radiation fields, the 



optically-thin parts of the stripped gas can be subject to the 
radiation drag and lose angular momentum. This might pro- 
vides a simple way to extract angular momentum from the 
interstellar medium. From a realistic point of view, however, 
the gas motion does not seem so simple. Recent high resolu- 
tion simulations on multiple supernova (SN) explosions in a 
star- forming galaxy (iMori. Umemura fc Fer rara 2004) have 
reveal that dense shells driven by SN hot bubbles undergo 
thermal instability and resultant cold clouds move together 
with ambient hot gas. Here, we consider two cases; one is 
a comoving case in which the stripping takes place solely 
by the radiation drag, and the other is a shear-fiow case in 
which the radiation drag is partially coupled with the ram- 
pressure stripping. 

The cloud temperature, mass, and radius are respec- 
tively assumed to be Tc = lO^K, Mc = 3.0 x W^Mq, and 
Rc = lOpc, while the temperature of hot ambient matter 
is Th = lO'^K. We postulate the peculiar velocity Vp of the 
cloud to be of the order of the virial velocity for a spheroidal 



galaxy with 
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As for the radiation fields, the radiation energy density is 
estimated by the assumed total luminosity L« and galaxy 
size to be 



E- 



^* - 3.45xl0"^erg cm"^ 



A-kRIc 
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3 X 1O12L0 / V Ikpc 



In this paper we assume the total luminosity to be L, 



3 METHOD OF CALCULATIONS 

We use the two-dimensional generalized curvilinear coor- 
dinates as shown in Fig. The number of grid points is 
101 X 101. Such coordinates allow us to take accurate bound- 
ary conditions for a spherical gas cloud and also to treat ap- 
propriately the flow of stripped gas on the right-hand side. 
The scale of computational domain is about 50pc. The gas 
cloud moves at a velocity Vp along the x-axis in the radia- 
tion fields. Here, to prevent the cloud from going out of the 
boundary with Vp, we make the Galilean transformation by 
Vp so that the cloud should stay at the origin if neither drag 
force nor ram pressure is exerted. In these coodinates, the 
surface of the cold gas is represented by a dash-dotted line 
in Fig. The other regions are filled with the hot gas. 

The two-dimensional Euler equations including the rel- 
ativistic radiation drag term can be written in the conser- 
vation form as 



dQ_ d£ dG 

dt dx dy 



H ^0, 



(3) 



where Q, F, G and H are respectively the following vectors. 
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G = 



H : 



pv 
puv 
pv^ + p 
(e + p)v 



PXE (u-vp) /c 
PXEv/c 




(4) 



where p is the density of the gas, u and v are respectively 
the X- and y- components of velocity, c is the speed of light, 
and X is the mass extinction coefficient due to dust opacity. 
In this work, we set x ~ lOOcm^g"^ by assumig the Solar 
abundance and 0.1/^m dust grains (e.g. see Umemura et al. 
1998). The pressure p and the total energy per unit volume 
e are related by the equation of state. 



(5) 



where 7 is the ratio of specific heats of the gas, which is 
chosen as 7 = 5/3. The radiation drag effect is taken into 
account by adding the term of H. 

To discretise the governing equation, eq @, the fi- 
nite volume-approach is used, where the AUSM-DV scheme 

is employed to ob tain the numerical convective fiux 

fWad a & Lioul ll994 . The time integration of a set of equa- 
tions is performed implicitly by a matrix free LU-SGS 
method (|Nakahashi et al.lll99<ill . 

The optical depth at each point is evaluated by inte- 
grating the opacity along with a line drawn from the outer 
boundary to the coordinate origin. The distribution of the 
optical depth in the cloud at initial time is shown in Fig. 
121 The optical depth at the center of cloud reaches above 
15 and the thickness of the optically-thin surface layer is 
~ 10% of the cloud radius. We involve the radiation drag 
term if the optical depth is lower than unity, because no drag 
force works in optically-thick medium iTsuribe fc Umemural 
[l993). 

In order to demonstrate the validity of the present 
method for the optical depth estimation, we show the opti- 
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Figure 2. The distribution of the optical depth in the vicinity 
of the cold gas at initial time. The optical depth at each point 
is evaluated by integrating the opacity along with a line drawn 
from the outer boundary to the coordinate origin. The surface of 
the cold gas is represented by a dash-dotted line. 



cal depth distributions on several snapshots of the dynami- 
cal evolution in Fig. |3] The dot-dashed semicircle shows the 
initial position of the cloud. The light gray or charcoal gray 
regions show the domain where the optical depth is from 
0.1 to 1 or above 1, respectively. The optically-thin regions 
are properly detected with tracing the dynamical evolution. 
The gas in charcoal gray regions is impervious to the ra- 
diation drag. However, there is limitation for applying this 
method. If the whole region of the cold gas deviates from 
the coordinate origin, the estimate of optical depth would 
not be correct. Threfore, at this time, the estimate of the 
radiation drag is not always accurate. In Section 2] we will 
be back to this point. 
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Figure 3. The optical deptli distributions on several snapsiiots 
of the dynamical evolution. The dot-dashed semicircle shows the 
initial position of the cloud. The light gray or charcoal gray re- 
gions show the domain where the optical depth is from 0.1 to 1 
or above 1, respectively. The gas in light gray is subject to the 
radiation drag, while that in charcoal gray regions is impervious 
to the radiation drag. 



4 RESULTS 

4.1 Comoving Case 

The stripping purely by the radiation drag is demonstrated 
in a comoving case, where no relative velocity is present be- 
tween the gas cloud and the ambient matter. The time varia- 
tion of density distributions is depicted by four snapshots in 
Fig. |11 and that of velocity fields is shown only for optically- 
thin cold gas in Fig. |3 where the density contours are also 



displayed. We can see in Fig.^Jthat the surface layers of the 
gas cloud are stripped by radiation drag and drifted to the 
right-hand side. Since we are observing the phenomena in 
the coordinates moving to the left at the velocity Vp, the gas 
drifted to the right is actually left behind by the radiation 
drag. It is noted that this stripping is obviously different 
from the ram-pressure stripping in that it occurs also be- 
hind the cloud in the downstream. The left-hand side of the 
optically-thin surface layer pushes the whole region to the 
right by the radiation drag. 

The velocity fields in Fig. I^show that the optically-thin 
surface layers move round the surface of the cloud. A part 
of the gas reaching the a;-axis is accelerated toward the right 
direction along the x-axis. This means that the stripped gas 
slows down by the radiation drag. 

In the later phase {t > lO^yr), the stripped gas is trailed 
by the radiation drag, because the cold gas is confined by 
the ambient hot gas. When the stripped gas reaches the x- 
axis, it splits into the left and right under the infiuence of 
external pressure. The radiation drag is not exerted on the 
gas splitting into the left owing to the large optical depth 
and therefore left there in the present Galilean transformed 
coordinates, while the radiation drag works to decelerate 
the optically-thin gas splitting to the right and eventually 
produces the gas flows to the right. As mentioned before, 
when the whole region of the cloud deviates from the origin, 
the present estimate of the optical depth in the left-hand 
side of the cloud becomes inaccurate. However, as far as the 
gas drifted to the right is concerned, the optical depth is 
still properly estimated. Hence, the trailing flow is thought 
to be real. 

4.2 Shear-Flow Case 

Here, we show the results for a shear-flow case in which the 
velocity difference between the gas cloud and the ambient 
matter is assumed to be O.lUp, corresponding to 0.13 in Mach 
number. The numerical results are shown in Fig. |S|and Fig. 
Q respectively for the density distributions and the velocity 
flelds for the optically-thin cold gas. 

The snapshots of the density contours in Fig. |S| show 
that the deformation of the gas cloud grows with time. In 
this case, the ram-pressure stripping occurs simultaneously 
with the stripping driven by the radiation drag. However, as 
shown in Fig. |S] (d) , the gas drifted to the right eventually 
flows along with the a;-axis. This flow pattern is similar to 
that in the case without the ram-pressure induced deforma- 
tion. 

The snapshots of the velocity fields show that the right- 
ward velocities of the optically-thin surface layers increase 
with the deformation. Our numerical method allows us to 
trace the optically-thin surface even in such strong deforma- 
tion of the cold gas. The optically-thin gas is subject to the 
radiation drag and makes a trailing fiow, similar to the case 
without the ram pressure. 

4.3 Momentum Loss 

The radiation drag works on both radial and azimuthal com- 
ponents of gas velocity, in proportion to each component. In 
a pure rotating system, the reduction only for azimuthal ve- 
locity leads to the angular momentum loss. But, in a chaotic 
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Figure 4. The time variation of density distributions in a co- 
moving case, (a), (b), (c) and (d) show these at time, 0.100 X 10^, 
0.388 X 10^, 0.504 x 10**, and 1.252 x lO^yr, respectively. The den- 
sity range is from 3.5 X 10~^^ to 5.3 X lO^^^g cm"''. Each density 
range is divided by 21 hnes with an equal increment. 



Figure 5. The time variation of vector fields for the optically-thin 
cold gas in a comoving case, (a), (b), (c) and (d) show the flow 
patterns at the same stages as Fig. |2 respectively. The density 
contours are also displayed. The density range and the number of 
contour lines are also the same as Fig. ^ 



system like a bulge under starbursts, the reduction of ran- 
dom velocities means, on average, the decrease of an az- 
imuthal component of velocity and therefore angular mo- 
mentum. In other words, the gas even in non-radial mo- 
tion can fall onto the center if the momentum of the gas is 
strongly reduced. Hence, here we consider the momentum 
loss in the present simulations. The momentum of optically- 
thin cold gas in the a;-direction, Qx, is defined as 



dV. 



(6) 



In Fig. |H1 the time variation of momentum is shown in units 
of the initial momentum. The solid line is the evolution in a 
comoving case, while the broken line is in a shear-flow case. 
A dotted line shows the level of 1/e reduction. 

From this result, we can evaluate the momentum loss 
time-scale, which is defined by the 1/e reduction of momen- 
tum. The time-scale for a comoving case \st = 1.04 x 10*yr, 
and that for a shear-flow case ist = 0.8 x lO^yr. In the shear- 
flow case, the surface gas of the cloud is stripped also by 
the ram pressure and is rarefied. Resultantly, the fraction of 
optically-thin gas increases compared to the pure radiation- 
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Figure 6. The time variation of density distributions in a shear- 
flow case, (a), (b), (c) and (d) show these at time, 0.100 X 10*, 
0.380 X 10*, 0.508 x 10*, and 0.862 x 10*yr, respectively. The 
density range and the number of contour lines are the same as 

Fig.ia 



Figure 7. The time variation of vector fields for the optically-thin 
cold gas in a shear-flow case, (a), (b), (c) and (d) show the flow 
patterns at the same stages as Fig. |5| respectively. The density 
contours are also displayed. The density range and the number of 
contour lines are the same as Fig.Bl 



drag case, so that the momentum is extracted more effi- 
ciently than the comoving case. 

The time-scale o f angular momentum extraction i s ana- 
lytically estimated bv lUmemura. Fukue fc Mineshig3 lll998l) 
to give 



2.4 X 10 yr 



fds 



3 x 1012L 



R„ 



1 kpc 



10- 



Ps 



°d 

0.1/im J \g cm 



(7) 



where /^^g, a^, and ps are the dust-to-gas mass ratio, the 
grain radius, and the density of solid material within the 
grain, respectively. It is found that the numerically obtained 
time-scale is in a good agreement with the analytic estimate, 
although the numerical one is slightly shorter. 



5 CONCLUSIONS AND DISCUSSION 

Wo have performed the radiation-hydrodynamic simulation 
for a gas cloud moving at a constant velocity in uniform ra- 
diation fields, with focusing attention on the radiation drag 
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Figure 8. The time variation of momentum in units of the initial 
momentum. The soHd Hne is the evolution in a comoving case, 
while the broken line is in a shear-flow case. A dotted line shows 
the level of 1/e reduction. 



as well as the coupling effect with the ram pressure. The 
obtained results are summarized as follows: 

(1) In a comoving case, the stripping occurs purely through 
the radiation drag. All the optically-thin surface layers are 
stripped from an optically-thick gas cloud, in contrast to the 
pure ram-pressure stripping. Finally, a trailing flow behind 
the cloud forms. 

(2) In a shear-flow case, the radiation drag-induced strip- 
ping occurs simultaneously with the ram-pressure stripping. 
But, the eventual flow pattern for optically-thin gas is sim- 
ilar to the pure radiation drag case. 

(3) The momentum loss time-scale for the stripped gas is 
on the order of 10* years under intensive radiation fields 
expected in an early evolutionary phase of a massive bulge. 
The time-scale is shorter for a shear-fiow case with the aid 
of the ram pressure. The present time-scale is basically in 
a good agreement with an analytic estimation for angular 
momentum extraction time-scale. 

It is demonstrated that the radiation drag by bulge 
stars is an effective mechanism to extract angular momen- 
tum from interstellar cloud in a galactic bulge system. This 
mechanism allows the mass accretion onto the galactic cen- 
ter, possibl y leading to the formation of a supermassive 
black hole fUmemura '2001'; 'Kawakatu & Umemura '2002'; 
teawa katu. Umemura. & Mori 2003; Kawakatu & Umemura 

In the present work, we have ignored the effect of self- 
gravity. The self-gravity is likely to be important for a cloud 
with Mc = 3.0 X lO^M© assumed in the present paper. Thus, 
an assumption of a uniform cloud may not be realistic, and 
instead the density distribution of a cold gas in gravitational 
equilibrium should be considered. In this case, the outside of 
a gas cloud could be rarefied and the domain of the optically- 
thin layers increases. Consequently, the efficiency of the ra- 



diation drag might be enhanced. To investigate the details 
of the stripping and momentum loss in such a case will be 
significant in the future work. 
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